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Abstract 



A nonlinear two-dimensional system is studied by making use of both the La- 
grangian and the Hamiltonian formalisms. The present model is obtained as a two- 

^ ■ dimensional version of a one-dimensional oscillator previously studied at the classical 

^ . and also at the quantum level. First, it is proved that it is a super-integrable system, 

H I and then the nonlinear equations are solved and the solutions are explicitly obtained. 

All the bounded motions are quasiperiodic oscillations and the unbounded (scattering) 
motions are represented by hyperbolic functions. In the second part the system is 
generalized to the case of n degrees of freedom. Finally, the relation of this nonlinear 
system with the harmonic oscillator on spaces of constant curvature, two-dimensional 
sphere 5^ and hyperbolic plane H^, is discussed. 
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1 Introduction and main results 

Mathews and Lakshmanan studied in 1974 [1],[2], the equation 

(1 + \x^)x- {\x)x'^ + a'^x = 0, A>0, (1) 

as an example of a non-hnear oscillator (notice a^ was written just as a in the original 
paper). In fact they considered (1) as a particular case of the differential equation 

y" + f{y)y" + 9{y) = o, 

that can be solved by using a two steps procedure: (i) first a reduction of order can be 
fulfilled by the change y' = p, y" = pp', p' = dp/dy; (ii) then the corresponding first order 
equation can be solved by using fi{y) = exp{2 J f (y) dy} , as an integrating factor. In this 
case we have f{x) = —\x/{l + \x^) and hence fj,{x) = 1/(1 + \x^); the general solution 
takes the form 

X = Asm{ut + 0) , 

with the following additional restriction linking frequency and amplitude 

„2 



UJ 
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1 + AA2 



That is, the equation (1) represents a non-linear oscillator with periodic solutions that they 
qualify as having a "simple harmonic form". The authors also proved that (1) is obtainable 
from the Lagrangian 



2 Vi + Aa;2 



which they considered as the one-dimensional analogue of the Lagrangian density 

L = l ( Y7^) id,<P 9^^ - ^" <P") , (3) 

appearing in some models of quantum field theory [3], [4]. 

The equation (1) is therefore an interesting example of a system with nonlinear oscillations 
with a frequency (or period) showing amplitude dependence. As a quantum system, the one- 
dimensional Schroedinger equation involving the potential x'^/{l + gx'^) was considered in [5] 
as an example of anharmonic oscillator, and later on studied in [6]-[10]; the three-dimensional 
quantum problem was considered in [11], [12] (e.g., in [11] the Bohr-Sommerfeld quantization 
procedure was applied in relation with some previous studies [13], [14], in nonpolynomial 
quantum mechanical models). We observe that this system can also be considered as an 
oscillator with a position-dependent effective mass (see [15] and references therein). 

The main objective of this article is to develop a deeper analysis of the equation (1) and 
the Lagrangian (2), first proving that this particular A-dependent nonlinear system can be 
generalized to the two-dimensional case, and even to the n-dimensional case, and second. 



pointing towards a simple geometric interpretation of the system so obtained. In more detail, 
the plan of the article is as follows: Sec. 2 is devoted to the properties of the A-dependent 
kinetic part of the Lagrangian and the A-dependent two-dimensional free motion. In Sec. 
3, that must be considered as the central part of this article, we study the A-dependent 
two-dimensional oscillator; we have divided this section in three parts: in the first part we 
discuss the existence of Noether symmetries for A-dependent Lagrangians of the form 

L{x, y, v^, Vy] A) = - {y~^) [^x + ^J + ^ i^'Vy - y^xf ] - V{x, y; A) . 
In the second part we discuss the properties of the A-oscillator described by the Lagrangian 



first proving that such system is superintegrable and then solving the equations of motion 
using the Lagrangian formalism as an approach; the third part deals with the Hamiltonian 
formalism and the A-dependent Hamilton- Jacobi equation which is shown to be separable 
in three different coordinate systems. In Sec. 4 we study the A-dependent n-dimensional 
nonlinear oscillator described by the Lagrangian 






X? 



and we solve the associated equations obtaining different types of solutions depending of 
the values of A. Next we consider the Hamiltonian approach, we prove that it is a super- 
integrable system and we obtain different families of constants of motion; the existence of 
several different sets of n commuting integrals is also discussed. In Sec. 5 we start with a 
discussion of Lie algebra structure of the symmetries of the n = 2 non-linear oscillator; in fact 
we see that they span a three-dimensional real Lie algebra isomorphic to 5*0(3, IR), 5*0(2, 1) 
or the Euclidean group in two dimensions, depending on the sign of the parameter A. Then 
we present a geometric approach which explains the surprising properties this A-dependent 
system has and relates it to the harmonic oscillator in spaces of constant curvature studied 
in Ref. [16]-[18]. Finally, in Sec. 6 we make some final comments. 



2 A-dependent "Free Particle" 

We will make use of the Lagrangian formalism as an approach. That is, we will look, in 
the first place, for a Lagrangian function L{X) with appropriate properties and then, we will 
turn our attention to the corresponding nonlinear equations arising from L{X). 

It is clear that the equation (1) represents a non-linear version of a linear equation with a 
non-linearity introduced by the coefficient A; but the important point is that, in Lagrangian 
terms, this coefficient A modifies not only the quadratic potential V = (1/2) x"^ of the 



harmonic oscillator but also the kinetic term T = (1/2) v^. Therefore, this particular system 
is not directly related with Henon-Heiles or any other similar non-linear system [19]- [22] 
where the nonlinearity is introduced by just adding a new term of higher order to the 
original potential. Our first aim is to extend to n = 2 dimensions this system in such a way 
that its distinguishing properties are maintained. So in order to construct the appropriate 
two-dimensional Lagrangian, we may split the problem in two: the problem of the kinetic 
term and the problem of the potential. 

Next we begin with that of the kinetic term. 

Before giving the expression for the two-dimensional kinetic term, let us consider some 
properties of the one- dimensional free-particle motion characterized by the following La- 
grangian 

2 Vl + Aa;2 
and the following equation 



L(x,t;.;A)=Ti(A) = -(^^^ , (4) 



(1 + Aa;2) X - (Ax) i;2 = , A>0. (5) 

Two important properties are: (i) The function Ti(A) is invariant under the action of the 
vector field X^ = Xx{\) given by 

d 



in the sense that we have 



X,(A) = Vl + Ax2 ^^ , 



X*(A)(Ti(A))=0, 



where X*(A) denotes the natural lift to the phase space IRxlR (tangent bundle in differential 
geometric terms) of the vector field ^^^(A), 

X*(A) = Vl + Xx^ — • ' "" ^ 



dx VI + Aa;2 ^ dv^ ' 

(ii) The general solution of the equation of motion, that can be directly obtained from the 
conservation of the energy E, is given by 

X = [-^) sinh (V2A^ (t + 0)) • (6) 

It is clear that this expression satisfies correctly the linear limit for A = 0. 

The first problem of the transition from 1-d to 2-d configuration space is the construction of 
the new two-dimensional T2(A). Many different A-dependent functions will have the same A = 
limit, so we must require that the new function T2(A) must satisfy certain properties. On 
the one hand we think that a natural requirement for T2(A) is to satisfy the two-dimensional 
versions of the previous points (i) and (ii). On the other hand the potential V^(A) for the 
two-dimensional oscillator (to be studied in the next section) must be a A-dependent central 



potential such that the angular momentum be preserved; but, according to the Noether 
theorem, exact symmetries of the potential lead to constants of motion only if they also are 
symmetries of the "kinetic energy". Hence, a necessary condition must be that T2(A) be also 
preserved with the same symmetry. 

We will consider as the starting point for our approach the following three requirements. 

1. The kinetic term T2(A) must be a quadratic function of the velocities that will remain 
invariant under rotations in the H^ plane. This means that it must depend of the 
coordinates x, y, by means of r'^ = x"^ + y'^, and of the velocities Vx, Vy, by means of 
vl + Vy, {xvy -yvxf, {xvx + yvyf and {xv^ + yvy){xvy -yv^). 

2. T2(A) should be invariant under (the lifts of) the two vector fields Xi{X) and X2(A) 
given by 

d 



Xi(A) = /TTaT^^, 

X2(A) = Vl + Ar2 I-, 

oy 

that represent the natural extension to IR^ of the vector field X^(A) associated to Ti(A) 
in the one-dimensional case. Notice that the previous point (1) implies that if T'2(A) is 
invariant under (the lift of) Xi{X) then so will be under (the lift of) X2(A). 

3. It must lead to a solution {x(t),y(t)) in IR^ for the two-dimensional free-particle mo- 
tion similar to the hyperbolic function x(t) in IR given by (6), for the one-dimensional 
free-particle motion. 

The most direct and simplest generalization of the one-dimensional kinetic term Ti(A) is 
given by 

T(x, ,,..,.,; A) = (i)(^^^^±^), 

but this function does not satisfy the point (2). So let us try a more general expression given 
by 

,1 W{Vx,Vy) 



T{x,y,Vx,Vy;X) 



■2' 1 + Ar2 ' 



where W = W{vx,Vy) denotes 



W = ci{vl + vl) + C2{xvy - yvxf + cz{xvx + yvyf + ca^xv^ + yvy){xvy - yv^) , 
then the point (2) is satisfied only if 

C2 = A Ci , C3 = , C4 = 0. 
Hence, we will choose the following two-dimensional kinetic function 

^2W = (^)(Y^77i)[^'+^S + ^(^^^-^^-)1' r' = x' + y\ (7) 



as the appropriate one for the two-dimensional A-dependent dynamics. Notice that this 
means that the A-dependence is introduced in two different ways: the first one is the global 
factor f/(l + Ar^) that is the most direct n = 2 extension of the one-dimensional factor 
1/(1 + Xx'^) in (2); the other A-term is not so simple and it represents a two-dimensional 
contribution that was not present (in fact, it can not be defined) in the one-dimensional 
case. Although one could guess that this additional term can introduce difficulties, we will 
see that it really simplifies most of properties, mainly all those related with symmetries (the 
geometric aspects will be discussed in Sec. 5). We also point out that we admit A can take 
both positive and negative values. It is clear that for A < 0, A = — |A|, the function (and the 
associated dynamics) will have a singularity at 1 — |A| r^ = 0; because of this we will restrict 
the study of the dynamics to the interior of the circle x'^ + y'^ < 1/|A| that is the region in 
which T2(A) is positive definite. 

It is known that a symmetric bilinear form in the velocities {vx,Vy) can be considered as 
associated to a two-dimensional metric ds"^ in H^. In this particular case, the function T2(A) 
considered as a bilinear form determines the following A-dependent metric 

ds^{X) = ( ^—-) \{l + Xy^)dx^ + {l + Xx^)dy^-2Xxydxdy] . (8) 

^1 + Ar"'^ L -I 

This relation between kinetic term and metric implies that the Killing vectors of the A-metric 
(is^(A) coincide with the exact Noether symmetries of the A-dependent free motion, that is, 
of the dynamics determined by assuming the kinetic term as Lagrangian, L{X) = T2(A). 

T2(A) remains invariant under the actions of the lifts of the vector fields Xi(A), X2(A), and 
Xj, given by 

Xi(A) = Vl + Ar2 ^, 

X2(A) = /TTaT^^, 

d d 

oy ox 

in the sense that, if we denote by X*, r = 1,2, J, the natural lift to the tangent bundle 
(phase space IR^xlR^) of the vector field X^, 

d , / xvr + yv,, \ d 






dy ^ Vl + Ar2 ^ dvy 

then the Lie derivatives of T2(A) with respect to X*(A) vanish, that is 

X*(A)(T2(A))=0, X*(T2(A))=0, r = l,2. 



We close this section by solving the two-dimensional free-particle motion determined by 
two-dimensional kinetic function (7). 
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The Euler equations arising from L{X) = T2(A), and representing the n = 2 generahzation 

of (5), are 



(1 + Ar^) X — A [i^ + y^ + A {xy — yxY] x = , 
{1 + \r^) y - \[x^ + y^ + \ {xy - yxf ] y = . 



(9) 



Of course they are much more difficuh of solving than the single equation (5); but assuming 
hyperbolic/trigonometric expressions for the two functions, x{t) and |/(t), then we obtain 
that the solutions of (9) are given by 

(_) sinh {Ct + 0i) , y= (-=) sinh {Ct + 02) 

A R 

' sin {Ct + 0i) , y= (^=) sin {Ct + ^2) 



X 



X 



A>0, 
A<0, 



(10) 



with the only restriction 



1, A>0, 
A^ + B^- A^B^ sin'^ id), -(])',) = 1, A<0. 



A' + B' + A'B' smh\(j)i - (j)2 

H - 02) 



After some calculus we arrive, for A > to 
hence the solutions {x{t),y{t)) reduce to 



\P^ = C\l-A^] 



C' = 2XE. 



X 



y 



j p2 _ \ J2 

'^^^ smh{V2XE {t + <!>,: 



IPj -\.P 
2XE 



sinh(V2AE(t + 02) 



:iii 



for A > 0, where J denotes the angular momentum xvy — yvx- For the case A < similar 
reasoning leads to 



X 



y 



Pf+ A|J2 


\ 2A|E 


Pi+ A|J2 


\ 2\\E 



sin A/2|A|^(t + 0i 



sin^/2|A|^(t + 02) . 



(12) 



These results generalize the hyperbolic solution x{t) in IR given by (6) and satisfy correctly 
the linear limit for A = 0. 



3 A-dependent n = 2 quasi-Harmonic Oscillator 



In this section we will study and solve the appropriate n 
equation (1) and the A-dependent Lagrangian (2). 



2 versions of the A-dependent 
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3.1 Noether symmetries of A-dependent potentials 

A general standard A-dependent Lagrangian (kinetic term minus a potential) will have the 
following form 

L{x, y, v^, Vy] A) = T2(A) - V{x, y; A) 

in such a way that for A = we recover the non-deformed linear system. 

It is known that if a potential V{x,y), defined the Euclidean plane, is invariant under ei- 
ther translations or rotations then it admits a Noether integral of first order in the velocities. 
Now, in this A-dependent case, we have also a rather similar situation, but as the Lagrangian 
system is A-dependent so are the transformations. The infinitesimal transformations gener- 
ated by Xi{\) are 



x' = X + e 6x , y' = y , Sx = vT+A 



r^ 



and can be interpreted as a A-dependent version of the translations along the x-axis. Similarly 
the generator of the one-parameter group of A-dependent translations along the y-axis 



X' 



' = X , y' = y + e 6y , 6y = Vl + Ar^ , 



is the vector field X2(A) (the generator of rotations remains A-independent). If we denote 
by 6l the Cartan semibasic one-form 

= ( . 2 ) v^dx + Vydy + \{xVy-yv^){xdy-ydx) , 

then we have the following 

1. If the potential V^(A) does not depend on the variable x then the Lagrangian L{X) 
is invariant under the transformations generated by Xi(A); in this case the function 
Pi (A) given by 



Vl + Ar2 

is a constant of motion. Notice that in this case the coordinate x is not cyclic since it 
is always present in kinetic term T2(A). 

2. If the potential V^(A) is independent of the variable y then the Lagrangian L{X) is 
invariant under the transformations generated by X2(A); in this case the function 
P2(A) given by 

P,(A) = ,:(.Y^(A))9, = J^=^ 

^ ^ V 1 + Ar^ 

is a constant of motion. This situation is similar to the previous one; that is, the 
coordinate y is not in the potential but, for At^O, it appears in the kinetic term. 



3. If V^(A) is a central potential, then 

J = i{Xj) 9l = xvy - yvx 
is a constant of motion. Notice that both the vector field Xj and J are A-independent. 

In these three very particular cases, the corresponding system becomes integrable with a 
second integral, -Pi (A), -P2(A), or J, arising from an exact Noether symmetry. 

3.2 Lagrangian approach 

Let us consider the following A-dependent Lagrangian 

L = T,{X)-V2{r;X), ^^^(r; A) = ^ (^^^) , r' = x' + y\ (13) 

where V2(r; A) is the direct extension to n = 2 of the n = 1 potential in (2). The dynamics 
is given by the following A-dependent vector field 

d d d d 

^\ = Vx^ + Vy — + F^{x, y, v^, Vy] A) h Fy{x, y, v^, Vy\ A) — 

where the two functions F^ and Fy are given by 



y 



in such a way that for A = we recover the dynamics of the standard 2-d harmonic oscillator 

^ d d / 2 \ "9 / 2 \ "9 

We see that for A < the potential V2(r; A) is a well with a boundless wall at r^ = 1/|A|; 
therefore, all the trajectories will be bounded. For A > we have that V2(r; A) — » (l/2)(a;2/A) 
for r -^ cxd; so for small energies the trajectories will be bounded but for E{X) > {l/2){a'^/X) 
the trajectories will be unbounded (see Figures I and II; notice that we have plotted V{x, A) 
but the graph of V2(r; A) is just the same but with r>0). 

Our objective is to solve the A-dependent equations arising from (13) and prove that this 
two-dimensional motion is periodic in the bounded case as it was in the one-dimensional 
case. At this point we recall that a system is called super-integrable if it is integrable (in 
the sense of Liouville- Arnold) and, in addition, possesses more independent first integrals 
than degrees of freedom; in particular, if a system with n degrees of freedom possesses 
iV = 2n — 1 independent first integrals, then it is called maximally super-integrable. An 
important property is that the existence of periodic motions is a characteristic related with 
super-integrability; thus we may suspect that this system is super-integrable. This is actually 
the case, and the following proposition states the super-integrability of this A-deformed 
system and proves the existence of a complex factorization. 

9 



Proposition 1 Let Ki, K2, be the following two functions 

X 



■ V 1 + Ar^ 
K2 = P2i\) + ia' ^ 



Vl + Ar2^' 
Then the complex functions Kij defined as 

K,, = K,K*, t,j = l,2, 

are constants of motion. 

Proof: We begin our analysis by considering tlie action of tlie vector field Tx, which represents 
the time-derivative, on the two A-dependent functions -Pi (A) and -P2(A). They are given by 



dt '' ' Vi + Ar2A^i + ^^2y' rft "' ' Vi + Ar2A^i + A^2y 



In a similar way, the time- derivative of the two velocity-independent functions, x/v 1 + Ar^ 
and y/V 1 + Ar^, is given by 



dt Vl + Ar2^ Vi + Ar2V ^^ ^ ' rft Vl + Ar2^ Vi + Ar^ 

Thus, the time-evolution of the two functions, i^i and 7^2, becomes 



dt dt dt ^Vl + Ar2 



and a similar calculus leads to 



lA,.r.(A,)^(^^)AV 



Thus we obtain 

d 



Therefore the potential ^(A) is super-integrable with the following three integrals of motion 

/i(A) = \K^\\ /2(A) = I i^2 P , h = Im(i^i2) = a {xVy - yv,) . 

That is, the existence of an invariant second order tensor Kij, admitting a complex factor- 
ization [23, 24], is preserved by the nonlinearity introduced by A. 

The system of equations 

(1 + Ar^) a: — A [i^ + y^ + A {xy — yx^] x + a'^ x = , 

10 



{1 + \r^)y- \[x^ + y^ + \{xy-yxy]y + a^y = 0, (14) 

cannot be directly solved in a simple way as it was the one-dimensional equation. Never- 
theless we can solve these equations by assuming certain particular expressions (with some 
undetermined coefficients) for the two functions x{t) and y{t). 

(i) Bounded motions: Let us look for solutions with the following periodic form 

X = A sm{ujt + (pi) , y = B sm{ujt + (f)2) , (15) 

where A, B, 0i, 02, and uj are real parameters. Then the equations (14) reduce to 

Ai?sin(c^t + 0i) = 0, BRsm{ujt + (f)2) =0, 
where R is given by 

i? = ^2 _ ^2 _ ^ (^^2 ^ ^2^ ^2 _ ^2^2^2 ^2 gj^2 ^^^ ^ ^^^ = 0^ _ 02 • 

Therefore the functions (15) are in fact solutions of (14) but with u, that represents the 
angular frequency of the motion, A-related with the coefficient a of the potential (that 
represents the frequency of the A = linear oscillator) by 

a^ = Muj^, M = l + APe, Pe = A'^ + B^ + X{A^B^sm'^(j)u)- 

Notice that the coefficient M is positive even for A = — |A| < since in that case the 
amplitudes A and B must satisfy A'^ + B^ < 1/\X\. 

Once we know the solution of the dynamics, we can obtain the constants values of the 
three integrals of motion, Ji, I2, an J; they are given by 

h = (I + A52 sin=^0i2)AV, 
h = (l + AA2sin2 0i2)fiV, 
J = —UJ AB sin 012 . 

Using these expressions, we can obtain the values of the amplitudes A, B, as functions of 
the integrals of motion; we have 

A = {-)^h-XJ\ B = {-)^h-\J\ 

UJ ^ UJ ^ 

so that the total energy becomes 

E{\) = i^){A' + B')co' + {^){ABsm<Pi2yuj\ 

< 



2^^1 + \PJ 2 A 

Let us summarize. The four coefficients (A, P, 0i, 02) remain arbitrary, the trajectories 
are sine-like periodic motions having the same frequency but (possibly) differing in amplitude 
and in phase, and the trajectories are "ellipses" in (x, y) plane. The situation is very similar 
to the one of the linear oscillator, the main difference laying in the frequency uj that is 
given by uj"^ = a'^/M and it depends, therefore, on the position (initial data). We have two 
possibihties: 

11 



• If the parameter A is negative A < 0, then oj > a. 

• If the parameter A is positive A > 0, then u < a. 

The energy can take any value for A < 0, and it is always bounded by E^^x = (l/2)(«^/A) 
for A > 0, with the value of u going down when the energy E[X) approaches to this upper 
value. 

(ii) Unbounded motions: Let us analyze the solutions corresponding to A > 0, and 

E > Ea,\. 

If we assume the following expressions 

x = A sinh(fit + 0i) , y = B sinh((]t + ^s) , (16) 

for x{t) and y{t), then we obtain that they are solutions of (14) with the condition that a 
and VL must be A-related by 

The constant values of the three functions, Ji, I2, an J, are given by 

/i = (1 + A 52 sinh^ 012)^2^2^ 
I2 = (l + XA^ sinh=^ 012) B^n^ , 
J = —Q AB sinh 0i2 , 

and the coefficients A, B, can be rewritten as follows 

A = i^)^h-\.P, B = i^)^h-\.P, 
Finally, the total energy becomes 

E{X) = i^)[A^ + B')n' + i^)[ABsmh(t>ufn^ 



2 ^ ^XPh- 1^ 2 A 

Figure III shows the form of the potential for A > and A < 0. 

(iii) Limiting unbounded motions: Let us analyze the very particular case character- 
ized by A > 0, and E = E^^x with E^^x = (1/2) (aVA). 

If we assume the following expressions 

x = Ait + Bi, y = A2t + B2, (17) 

for the solutions, then we obtain that they are solutions of (14) with the following A- 
dependent restriction for the four coefficients Ai, A2, Bi, and B2, 

a^ = \PL, PL = Al + Al + \{A2Bi-AiB2f. 

12 



In this particular case the three functions, Ji, /2, an J, take the form 

h = Al + X{A2B^-A,B2f, 
h = Al + \{A2Bi-AiB2)\ 
J = A2B1-A1B2, 

and the coefficients Ai, A2, can be rewritten as follows 

A^ = ^h-X.P, A2 = ^l2-\-P. 
Finally, the total energy becomes 

E{X) = (i) (Aj + Al) + (^) {A2B, - A,B2f = ^ 



2' V ^ ' ^y ■ ^2' ' ' 2A 

Thus, the linear functions (17) appear, in this nonlinear system, as a border-line solution 
making separation between two different situations in the A > case: trigonometric periodic 
oscillations (15) for small energies and hyperbolic unbounded (scattering) evolutions (16) for 
high energies. 



3.3 Hamiltonian approach 

The Legendre transformation is given by 

_ (1 + \y'^)vx - Xxyvy _ (1 + \x'^)vy- Xxyv^ 

Px ~ 1 I \ 2 1 Py ~ 1 I \ 2 ' 

l + Ar"' 1 + Ar"^ 

so that the form of the angular momentum is preserved by the Legendre map, in the sense 
that we have xpy — ypx = xvy — yv^ (notice that this fact is consequence of the introduction 
of the term J in the definition of T2(A)), and the general expression for a A-dependent 
Hamiltonian becomes 



.1. 



'2' 



H{x,y,Px,Py;X) = {-) p^ + Py + X (xp^ + yPy) +i-)a V{x,y) (1 



1, 



•2' 



and hence the associated Hamilton- Jacobi equation takes the form 
^dS\2 fdS\2 ^ / dS dS\2 2t./ n 

This equation is not separable in (x, y) coordinates because of the A-dependent term; indeed 
the {x, y) coordinates are not even orthogonal. Nevertheless we will see that there exist 
three particular orthogonal coordinate systems, and three particular families of associated 
potentials, in which the Hamiltonian (18) admits Hamilton- Jacobi separability. To find 
them, we look for two one-parameter families of curves fi{x,y,Ci) = and f2{x,y,C2) = 0, 
(?( A) -orthogonal to the level set curves of x and y, we obtain respectively. 



x = Ciyl + Xy'^, and y = C2V 1 + Ax^ . 
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These two expressions suggest us to consider two particular systems, that we will denote by 
{zx,y) and {x,Zy), that can be seen as A-dependent deformations of the {x,y) coordinates; 
the third system is just the polar coordinate system. 

(i) In terms of the new coordinates 

V^ + ^y 

the Hamilton- Jacobi equation becomes 

(^ + ^-l){^y + (^ + ^yr {^y + a^l + ^y')V = 2(1 + Xy^)E 
so if the potential V{x,y) can be written on the form 

V=ifi(f£l + M/,fe) (19) 

1 + Ay^ 

then the equation becomes separable. The potential is therefore integrable with the following 
two quadratic integrals of motion 



^ ^ 1 + Xy^ 

fii) In terms of 



( \ y 

yX, Zy) , Zx 



V 1 + A X" 



which is the symmetric one of the previous change (i), the Hamilton- Jacobi equation becomes 

(1 + Ax2)2(|^)' + (1 + a4)(|^)' + a^ (1 + Xx')V = 2(1 + Xx')E, 
and, therefore, if the potential V{x,y) can be written on the form 

then the equation becomes separable. The potential is therefore integrable with the following 
two quadratic integrals of motion 

/i(A) = {l + Xr^)pl-XJ' + a'(w,{x)--^^^W^{zy)), 



1 + Aa;2 



/2(A) = {l + Xr')p'y + a'W2 



-yj- 



(iii) Finally we may use polar coordinates {r,(f)). Here the A-dependent Hamiltonian (18) 
is given by 

H{r, <P,pr, p^; A) = (^) [ (1 + A r')pl + ^] + (^) \/(r, 0) 
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so that the Hamilton- Jacobi equation is given by 



i+^^^)(f)'4(ir+"^''(^»-2^- 



dr ■' r 
Let us suppose that the potential V takes the form 

G(0) 



V = Fir) 



..2 ' 



(21) 



then the equation admits separability 

'dSx2 



r^(l + Ar2)(:^)Vr2(a2F(r)-2E)] + [(|^)' + a2G(0) 



0. 



The potential V is integrable with the following two quadratic integrals of motion 

l-r2 



(l + Ar2)p2+(^L)p2^«2 



/i(A) 

/2(A) = 4 + «2G(0). 



Fir) 



G(0) 



It is clear that f G = then V is a A-dependent central potential and the function I2 just 
becomes I2 = p^. 

In these three separable cases the potential V(A) is integrable with two quadratic constants 
of motion, /i(A) and /2(A), in such a way that 



^(A) = d)(/i(A) + /2(A)) 



That is, there exist three different ways in which the Hamiltonian FLiX) admits a decompo- 
sition as sum of two integrals; notice that in the linear case A = they reduce to only two 
since the two first cases, (i) and (ii), coincide. 

We recall that a potential V is called super-separable if it is separable in more than one 
system of coordinates [25]-[30]. The potential 



can be alternatively written as follows 



a 



V^2(A) = — 



2 Vi + A(a;2 + |/2) 



1 



2 Vi + Ay2yLi + Az^ 



a 



1 



2 Vi + Aa;2 



x^ + 



1 + Az2 



a 



2 Vi + Ar2 
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Therefore, it is super-separable since it is separable in three different systems of coordinates 
[zx,y), {x,Zy), and {r,(f)). Because of this the Hamiltonian 



H = (-) [pI + pJ + a {xp^ + ypyf 



a 



x^ + y^ 



2 Vi + A(x2 + |/2) 



(22) 



admits the following decomposition 

H = Hi + H2 — XH^ 
where the three partial functions i^i, H2, and H^, are given by 

1 



Hi 

H2 
H3 



2 

1 r 



X 



{l + Xr')pi + a ,^^,2 
{xpy - yp^f , 



each one of these three terms has a vanishing Poisson bracket with H for any value of the 
parameter A 

{H,Hi} = Q, {H,H2} = 0, {H,H,} = 0. 

So the total Hamiltonian can be written as a sum, not of two, but of three integrals of 
motion. The third one, which represents the contribution of the angular momentum J to 
H, has the parameter A as coefficient; so it vanishes in the linear limit A — * 0. 



4 n-dimensional Oscillator 

We have studied with a great detail the A-dependent nonlinear bi-dimensional oscillator; 
nevertheless, it is clear that this particular Lagrangian (Hamiltonian) system admits a direct 
extension to n dimensions. 



4.1 Lagrangian formalism 

The n-dimensional A-dependent Lagrangian is given by 



a 



Kj 



2 Vi + Ar^ 



(23) 



where we have made use of the notation 



7 Jb : • CXllLi '-^ %') "^ Z I i % 5 5 / X....»/i. 
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The A-dependent Euler-Lagrange vector field Tx takes the form 

d d 

Ta = Vk ^ h Fk{x, V, A) -— 

OXk OVk 

where the n functions F^ = Fk{x, v, \), k = 1, . . . ,n, are given by 



^' = -Hrfb) + ME^'? + ^ E^Krfb) ■ 



Same as in the n = 2 case, we can obtain the exphcit expressions for the solutions of the 
dynamics. We also have, in this n-dimensional case, a quasi-harmonic oscillatory motion for 
A < 0, and two different qualitative behaviours with a border-case when A > 0. 

(i) Periodic motions: We can assume that the solutions Xi = Xi{t) of the n equations 

Xi J^iyX^ X, A J , t 1, . . . , Tt, 

are periodic functions of the form 

Xi = Ai sm{uj t + (f)i), i = 1, . . . , n. 

Then we arrive at 

Ai R sin(a; t + (pi) = , i = 1, . . . ,n, 

where R is given by 



R = a^ - (l + XJ2A^) uj^ - X^ (j^AlA^sin^ (j)^. 



Jij 



ijj . 



*j 



Therefore, we have uj^ 7^ a^ but uJ^ = a'^/M with M = M(A) given by 
M(A) = 1 + A P^-^ , P^-^ = ^ A^ + A ( E A'4 sin^ 






and the energy E{X) is given by 



Eix) = ih {j: ^) ^' + d) J:{^^ A, sin <i>,,yu' 



a\ / P!;^^ 



2^h + APi")^' 

(ii) Unbounded motions: If we assume that the solutions Xi = Xi{t) are of the form 

Xi = Ai sinh(i7 1 + (pi) , i = 1, . . . ,n, 
then we arrive to Q^ = a'^/M with M(A) given by 

M(A) = -1 + A Pt^ , Pi"^ = E ^' + ^ ( E 4'^? sinh^ 0,, ) , 
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and the energy E[X) is given by 

E{X) = {b{J:^)^' + d)T.{^^AJsm<j>,,yn' 



«J 



I 

(iii) Finally, in the positive A > case, there exists a very particular border-case that 
makes separation between the behaviours (i) and (ii). It is characterized by the a value of 
the energy E given hy E = E^^x, and its time-evolution is represented by linear functions 

Xi = Ait -\- -t>i , i = 1, . . . , n, 

with the following A-dependent restriction for the 2n coefficients Aj, Bi, 

a' = \PL, pP = E ^' + ^ Y.(A^B, - A.Bif . 

Concerning the energy E{X) it takes the particular value 

E{\) = {\) E A' + (^) E(A^. - A3B^f = ^ . 

4.2 Hamiltonian formalisin 



The n-dimensional A-deformed Hamiltonian is given by 

2" Vl + Ar2 



^ = ^[Ep^ + A(E-^p01+^(t-^), r^ = T.-l (24) 



I 



It can be rewritten as a sum of A^ = (l/2)n(n + 1) quadratic terms as follows 



^ = d)(E4(A))-d)(E4 



where the functions /a;(A) and Jjj, which are given by 



xl 



4(A) = {l + Xr')pl + a\-^^), A;=l,...,n, 



+ 

*i XiPj XjPi , i, J i, . . . , vi, 

are constants of motion 

{H,lki\)} = 0, {H,J,^} = 0. 



Moreover, the A-dependent functions 

I,,{\) = {l + \r^)p,p, + a'{j^^), t,j = l,...,n, (25) 

are constants of motion as well. Of course, the existence of these three different families, 
/fe(A), Jij, and Iij{X), means that all these integrals cannot be independent since the max- 
imum number of (time-independent) functionally independent integrals is A^ = 2n — 1. In 
order to obtain a fundamental set of independent integrals we can choose the n functions 
/fe(A) and n — 1 of the angular momenta; an example is given by 

(/fc(A) , Ji,i+i) , k = l,...,n, i = l,...,n-l. 

In fact this situation is rather the same that one finds in the linear A = case but with 
two important distinctions: first that all these facts remain valid also for the unbounded (or 
scattering) motions, present when A > and second that the algebra of Poisson brackets 
seems to be quadratic. 

We close this section observing that, if we call super-separable to a system that admits 
Hamilton- Jacobi separation of variables (Schroedinger in the quantum case) in more than 
one coordinate system, then quadratic super- integrability (i.e., super-integrability with lin- 
ear or quadratic constants of motion) can be considered as a property arising from super- 
separability. In this A-dependent case, the three families of constants of motion are of such a 
class {Jij are linear and /^(A) and Iij{\) are quadratic), so we can conclude that the super- 
integrability of the Hamiltonian (24) arises from its multiple separability. In the linear A = 
case, since the Hamiltonian is directly separable in the n-dimensional cartesian system, the 
Hamiltonian is just the sum of the n partial one-dimensional energies; in the general non- 
linear At^O case, the functions /j(A), and Jij i,j = 1, . . . ,n, arise from separability in the 
n-dimensional versions of the two-dimensional coordinates {zx,y) and {x,Zy) studied in Sec. 
3 and in the n-dimensional spherical system. Finally, there exist many different sets of n 
conmmuting constants; as an example, in the n = 3 case we have the following three sets of 
involutive integrals 



/i , J2 — A Ji2 , ^3 — A ( J23 + J- 
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ih-\iJj, + Jl),h,Is-\Jl,), 

( Jl — A Jg]^ , J2 — A ( J^2 + -^23) ) -^3 ) ; 

as well as (/1 + /2 + /3, Jij, J^). It is clear from this example that in the general n-dimensional 
case there exist many more different ways of constructing involutive sets of n integrals. 



5 A Geometric Interpretation 

In this section we will discuss some additional properties that will prove to be related with a 
new geometric interpretation. In particular, we will see that this quasi-harmonic nonlinear 
oscillator turns out to be closely related with the harmonic oscillator on a space of constant 
curvature. 
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5.1 Some additional properties 

We begin by considering two remarkable properties; the first one is concerned with the 
symmetries of two-dimensional system and the second one with the Lagrangian of the one- 
dimensional oscillator. 

(i) The three A-dependent vector fields Xi, X2, Xj, obtained in Sec. 2 as the Killing vectors 
of the metric ds'^ or, equivalently, as the Noether symmetries of the Lagrangian L{X) = T2(A) 
of the A-dependent free particle, close the following Lie algebra, 

[Xi(A) ,X2(A)] = \Xj, [X^{\),XJ] = X,{\) , [X^W^Xj] = -Xi(A) . (26) 

This means that the 2-d configuration space has a three-dimensional symmetry Lie algebra, 
hence a maximal one, and this implies that the space should be of constant curvature. Indeed 
the Lie algebra (26) is isomorphic to the Euclidean algebra in the particular A = case, and 
to the Lie algebra of the isometries of the two-dimensional spherical (A < 0) and hyperbolic 
spaces (A > 0) in the general A 7^0 case. 

It seems, therefore, that there exists a certain relation between this A-dependent nonlinear 
oscillator and the properties of the two-dimensional spaces of constant curvature. At this 
point we recall that we have proved in Sec. 4 that this nonlinear system is well defined for 
any number n of degrees of freedom; therefore, it is natural to guess that if such relationship 
exists then it must be true, not only for n = 2, but for any dimension. 

Although the n = 1 case can be considered as a very special case, it seems convenient 
to go back to the one- dimensional oscillator and analyze again its properties but now in 
relation with the above point (i). We have obtained the following property that concerns 
the one-dimensional nonlinear system. 

(ii) Let us consider the change (x, w^;) —^ ilyVg) given by 



q = (—1=) sinh ^{VXx) , A > , 



then the Lagrangian 



becomes 



L{x, v^; A) = i (y^) i^l - «' ^') (27) 



a^ 



L{q, V,- A) = - vl - (^— j tanh2(v^g) . (28) 

In the negative case, A < 0, A = — |A|, the corresponding change is given by 

g=(-L)sin-i(V^x), A = -|A|, 
V 1^1 



and then we arrive at 

Tin 7).- X^l ^ — -J, 



L(g, V,- \) = ]-vl-[—-] t^n\J\\\ q) . (29) 
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Hence, we can remove the A parameter from the kinetic energy Ti(A), but the price for this 
simphfication is that the potential V drops its rational character and becomes a trigonometric 
or hyperbolic squared tangent function 

Viq;\) = ^(y)tanh2(VAg), for A>0, 
V{q;\) = -(^)tan2(y^g), for A<0. 

From a practical point of view, this change does not simplify really the problem, since the 
two new equations 

./\X\q + a^F{q;X) =0, 

with F = F{q, A) given by 

sinh(v^g) , , ^ sm{J\X\q) 

F = ; V , ^ > , and F= ^ ^ , A < , 

cosh-'(VAg) cos%^\\\q) 

are not easier of solving than the original one (1). Nevertheless, from a more qualitative 
viewpoint, the new aspect adopted by this one-dimensional Lagrangian is a very interesting 
fact since the form obtained for the new potential V{q, A) is closely related with the polar 
coordinate (p, 0) expression for the two-dimensional potential of the harmonic oscillator on 
two-dimensional spaces of constant curvature, previously studied (by two of the present 
authors) in Refs. [16]-[18]. 

It is clear from properties (i) and (ii) that we can conjecture the existence of a direct 
relationship between this particular nonlinear oscillator and the harmonic oscillator on spaces 
of constant curvature. Another remarkable property is that the set of integrals of motion 
(25) is analogous to the ones appearing in the n-dimensional version of the Smorodinsky- 
Winternitz system in curved spaces [31]. 

5.2 The harmonic oscillator on spaces of constant curvature 

Higgs [32] and Leemon [33] analyzed the characteristics of the two fundamental central 
potentials, Kepler problem and harmonic oscillator, on the A^-dimensional sphere. Since 
then a certain number of authors have studied this question from both the classical and 
the quantum points of view [34]-[40]. Next we recall some of the basic properties of the 
formalism studied in [16]-[18]. 

We begin for the following three basic ideas: 

• The harmonic oscillator is a system that is well defined in all the three two-dimensional 
spaces of constant curvature (sphere S'^, Euclidean plane E^, and hyperbolic plane H"^). 
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• A joint approach, where the usual Euchdean system and the two curved systems (de- 
fined on S*^ and H^) can be studied, all the three, at the same time, is possible using 
the curvature /t as a parameter. 

• Thus the spherical and hyperbolic oscillators can be considered as curvature deforma- 
tions of the well known "fiat" Euclidean oscillator which arises as a very particular 
case of the more general "curved" systems. 

This two last points mean that all the theory must be developed by making use of k- 
dependent mathematical expressions leading to general /t-dependent properties. The specific 
properties characterizing the harmonic oscillator on the sphere, on the Euclidean plane, or 
on the hyperbolic plane, are then obtained particularizing for k, > 0, k, = 0, or k, < 0. In 
order to present these expressions in a form which holds simultaneously for any value of k, 
the theory can be developed by making use of the following "tagged" trigonometric functions 

cosv^ft^a; if /t > 0, i^^^^\^^ if k > 0, 

Cn{x) = { 1 if K = 0, S«(a;) = I x if k = 0, (30) 

coshv^^a; if /t < 0, [ 7= sinhv/^x if k < 0, 

If we make use of this /t-dependent notation in the polar coordinates (p, 0) system, then the 
differential element of distance, on the three spaces (S*^, E , H'^), can be written as follows 



It reduces to 

dsl = dp^ + (sin^ p) dcp^ , ds^ = dp^ + p^ dcj? , ds^_^ = dp^ + (sinh^ p) dcp^ , 

in the three particular cases of the unit sphere, Euclidean plane, and 'unit' Lobachewski 
plane. Consequently, the Lagrangian for the geodesic (free) motion on the spaces (S*^, E^, H'^) 
is given by the kinetic term arising from the metric 



L{k) = TiK) = (-) [v^^ + sUpH 

and the Lagrangian for a general mechanical system (Riemmanian metric minus a potential) 
has the following form 



L{'^) = ih {vl + tipH)-Uip,<l>,'^) 



'2' 
It is clear that the well known expression for a natural Lagrangian on the Euclidean plane 

L = (^) (vl + / vl) - Vip, 0) , Vip, 0) = f/(p, 0, 0) , 

is obtained as the particular k, = case of L(/t). 
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Until now we have considered general aspects of the theory of Lagrangian systems on 
curved spaces. Now let us turn our attention to the spherical and hyperbolic harmonic 
oscillators; it is characterized by the following Lagrangian with curvature k 



{\){vl + ^l{p)vf)-{\)ujlTl{p). 



where the /t-dependent tangent Tk(p) is defined in the natural way Tk(p) = Sk(p)/Ck(p)- 
In this way, the harmonic oscillator on the unit sphere (Higgs oscillator), on the Euclidean 
plane, or on the unit Lobachewski plane, arise as the following three particular cases 

f/i(p) = (i)c.o^tanV, V{p) = U,{p) = {]^)ulp\ f/_i(p) = (i) Cq^ tanh^- 

The Euclidean oscillator V{p) = Uo{p) (parabolic potential without singularities) appears 
in this formalism as making a separation between two different situations. The spherical 
potential is represented by a well with singularities on the border (impenetrable walls at the 
equatorial circle p = 7r/2y^ if the potential center is placed at the poles), and the hyperbolic 
potential by a well with finite depth since for /t < 0, /t = — |/t|, we have Ui^{p)^{1/2)[u!q/\k\) 
when p -^ oo. Actually, the Scarf potential V{x) = 7/cos^(a;), which differs in a constant 
term from Ui{x), has been studied in solid state physics and has many interesting properties 
[41]- [43]. Figure IV plots U^ip), for the three particular cases of the unit sphere, Euclidean 
plane, and 'unit' Lobachewski plane; notice the great resemblance with Figure III. 



5.3 On the existence of a relation between two dynamics 

The abovementioned property (ii), for the n = 1 case, suggests that the Lagrangian (27) 
and the equation (1) is a nonlinear model for an harmonic oscillator in the circle S^ (A < 0) 
or in the hyperbolic line (A > 0). In a similar way, if we consider the two properties (i) and 
(ii) together, we arrive at the conclusion that this correspondence must also exist for the 
nonlinear n = 2 oscillator. Remark that in the n = 1 case, S^ and H^ must be understood 
as one-dimensional spaces obtained by endowing each single geodesic of S'^ or H^ with the 
induced metric. Motion on S^ and H^ will correspond to the J = radial motions on S'^ or 

Hence, the Lagrangian (27) for the equation (1) can be considered as a nonlinear model 
on the IR-line for the harmonic oscillator on the circle S^ and the hyperbolic line H^; the 
Lagrangian L = T2(A) — V2(A) of equation (13) as a nonlinear IR^-model for the harmonic os- 
cillator on the sphere S"^ and the hyperbolic plane H^; finally, the more general n-dimensional 
Hamiltonian given by equation (24) as a nonlinear lR"-model for the harmonic oscillator on 
the n-dimensional spaces S*" and iJ". So we arrive to the conclusion that an harmonic 
oscillator, that is a linear system, when is defined on a space of constant curvature turns 
out to be equivalent to a nonlinear oscillator on IR, IR^ or IR", with the nonlinear parameter 
A playing the role of the (negative of the) curvature k. Notice that, in dynamical terms, 
this equivalence is a relation of conjugacy, and also that the existence of this relation can be 
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considered as the origin of the quasi-harmonic behaviour obtained for the solutions of the 
nonhnear system. 

The existence of this relationship makes this nonlinear H^ (or IR") system even more 
interesting. From an abstract or qualitative viewpoint, the system on the sphere S'^ (or on 
the hyperbolic space H^) may be considered as the more fundamental one; nevertheless, we 
have been able to solve the equations and to obtain the explicit solution for the dynamics 
because we were working with the nonlinear A-dependent H^-model (or lR"-model). Hence, 
the nonlinear system studied in this article results to be much more appropriate for the 
explicit resolution of problems. 



6 Final Comments and Outlook 

We have proved that the Lagrangian (13) is the appropriate two-dimensional version of the 
Lagrangian (2), (27), and we have solved the nonlinear equations (14) in the two cases 
A < and A > 0. Moreover we have proved in Sec. 4 that this nonlinear system admits a 
n-dimensional version given by (23) and (24), and we have also solved the correspondent 
system of n equations. 

We think that all these results suggest the study of some related questions among them 
we point out the following: Firstly, the Poisson brackets of the integrals of motion for 
the 2-dimensional Hamiltonian (22) or the n-dimensional Hamiltonian (24) seems to close a 
quadratic Poisson algebra (see e.g. Ref. [44]); this is a very interesting problem that deserves 
to be studied. Secondly, the geometric relations obtained in the last Sec. 5, is a matter 
related with the existence of two different but conjugate dynamical systems. Conjugate 
systems are systems related by diffeomorphisms preserving the fundamental properties; so 
it is convenient to develop a deeper analysis of this particular relation of conjugacy existing 
between the nonlinear A-dependent oscillator and the harmonic oscillator on the sphere S"^ 
or the hyperbolic plane H^. Thirdly, it is clear that the calculus are easier of handle when 
working with the A-dependent formalism than when working on the constant curvature 
spaces; therefore, when looking for new dynamical results on the spaces such as S"^ or H^, 
a practical strategy will be to first consider the question by using Lagrangians such as (13) 
or (23). Fourthly, in Refs. [17] and [18] we have studied not only the central 1:1 oscillator 
given by the function Vk(p) = (1/2) cJq Tk(p) but also some others non-central non-isotropic 
oscillators Vk(p, 0); it will be convenient to consider these non-central problems also inside 
this A-dependent formalism. Finally, we mention the study of the quantized versions of all 
these nonlinear systems. Notice that, from the quantum point of view, as it is an oscillator 
with a position-dependent effective mass, there is a problem with the appropriate order for 
the factors in the kinetic term. We think that all these problems are examples of some open 
questions that must be investigated. 
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Figure Captions 

Figure I. Plot of V{\) = (1/2) (a^x^)/(l + Ax^), a = 1, A < 0, as a function of x, for 
A = —2 (upper curve), and A = — 1 (lower curve). 

Figure II. Plot of V{\) = (1/2) (a^a;^)/(l + Ax^), a = 1, A > 0, as a function of x, for 
A = 1 (upper curve), and A = 2 (lower curve). 

Figure III. Plot of ^(A) = (1/2) (aV)/(l + Ar2), a = 1, as a function of r, for A = -1 
(upper curve), A = (dashed line) and A = 1 (lower curve). 

Figure IV. Plot of Uf,{p) = (1/2) cj^ T^(p), cjq = 1, as a function of p, for k = -1 (lower 
curve), K = (dash line), and k, = 1 (upper curve). 
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Figure I. Plot of V{X) = (1/2) (aV)/(l + Xx^), a 
= —2 (upper curve), and A = — 1 (lower curve). 
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Figure II. Plot of V{X) = (1/2) (aV)/(l + Xx^ 
A = 1 (upper curve), and A = 2 (lower curve). 
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Figure III. Plot oiV2{X) = (1/2) (aV)/(l + Ar2), a = 1, as a function of r, for A 
(upper curve), A = (dashed line) and A = 1 (lower curve). 
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Figure IV. Plot of Uf,{p) = (1/2) lu^ tUp), ujq = 1, as a function of p, hi k = -I (lower 
curve), K = (dash line), and k, = 1 (upper curve). 
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